Quantum-classical transition of the escape rate of uniaxial antiferromagnetic particles 

in an arbitrarily directed field 



o 
o 

(N 
> 

o 



Bin Zhou 1,2,3 , Ruibao Tao 1 , and Shun-Qing Shcn 2 
1 Department of Physics, Fudan University, Shanghai 2004-33, China 
2 Department of Physics, The University of Hong Kong, Hong Kong, China 
A Department of Physics, Hubei University, Wuhan 430062, China 
(Dated: February 2, 2008) 

Quantum-classical escape rate transition has been studied for uniaxial antiferromagnetic particles 
with an arbitrarily directed magnetic field. In the case that the transverse and longitudinal fileds 
coexist, we calculate the phase boundary line between first- and second-order transitions, from 
which phase diagrams can be obtained. It is shown that the effects of the applied longitudinal 
magnetic field on quantum-classical transition vary greatly for different relative magnitudes of the 
non-compensation. 
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Escape from a stable or metastable state at high tem- 
peratures is governed by a classical thermal activation 
rate. At low temperatures close to zero, quantum tunnel- 
ing becomes relevant. When these two escape rates are 
equal there exists a crossover temperature To at which 
a transition between classical and quantum regimes oc- 
curs. The study of the quantum-classical transition is 
an interesting subject with a long history 1] One of the 
main issues in this subject is to determine whether the 
transition is first- or second-order. The transition was 
recognized as a smooth second-order one in the quan- 
tum mechanical models of Afflecky! and the cosmologi- 
cal models of Linde. 2j However, it was shownQ that the 
smooth transition is not generic. Chudnovsky has sug- 
gested that the order of transition is determined by the 
behavior of the Euclidean time oscillation period t(E), 
where E is the energy near the bottom of the Euclidean 
potential, which corresponds to the top of the potential 
barrier. |3| The non-monotonic behavior of the oscillation 
period as a function of energy, i.e., the existence of a 
minimum in the t ~ E curve, was proposed as a condi- 
tion for the first order transition in quantum-mechanical 
tunneling. [3| Later, a sufficient criterion for the first- 
order phase transition was obtained by carrying out the 
nonlinear perturbation near the sphaleron solution. £| 

Since the first- and second-order transitions between 
the quantum and classical behaviors of the escape 
rates in spin systems were introduced by Chudnovsky 
and Garaninjajgthe topic has attracted considerate 
attention. 0, UlallS 13 111 13 13 El Most theoretical 
studies have been focused on the ferromagnetic particles. 
However, most ferromagnetic systems are actually ferri- 
magnetic particles. For instance, both Mn^Ac and Fe§ 
are characterized by a large spin ground state which orig- 
inates from incomplete compensation of antiferromag- 
netically coupled spins. |16( The strong exchange interac- 
tion should be taken into consideration. In Rcf. [13 Kim 
treated the phase transition in ferrimagnetic or antifer- 
romagnetic particles for two general forms of the mag- 
netic anisotropy energy. Very recently, the quantum- 



classical transition in antiferromagnetic particles with 
biaxial symmetry in the presence of an applied mag- 
netic field along the medium axis or along the easy axis 
was investigated. |18lll9| Note that recent work of Chud- 
novsky and Garanin postulates dislocations as the main 
source of spin tunneling in Mni2 crystals. [2(j Their the- 
ory shows that when the external magnetic field is ap- 
plied along the c axis of the crystal, local rotations of the 
magnetic anisotropy axis due to dislocations result in the 
effective local transverse magnetic field. Experimental 
evidence of the effects of dislocations on tunneling has 
been also reported. |2lJ, l22l l2.il 1 12 Tl ] Therefore in the study 
of the quantum-classical transition of Mni2 , the case of 
coexistence of the transverse and longitudinal magnetic 
field is worth investigating. Considering the molecular 
cluster Mni2 actually is ferrimagnetic, the exchange in- 
teraction should be also taken into account. In this paper 
we aim to investigate the quantum-classical transition of 
the escape rate of uniaxial antiferromagnetic particles in 
an arbitrarily directed field, i.e., the coexistence of the 
transverse and longitudinal magnetic field. It is shown 
that the effects of the applied longitudinal magnetic field 
on the quantum-classical transition vary greatly for dif- 
ferent relative magnitudes of the non-compensation. 

We consider a small uniaxial antiferromagnetic particle 
with two magnetic sublattices whose magnetizations, mi 
and m2, are coupled by the strong exchange interaction 
mi-ni2/xj_, where x± is the perpendicular susceptibility. 
The system of interest has a non-compensation of sublat- 
tice with m (= m% — rri2 > 0), and easy-axis anisotropy 
along the z-axis. In the presence of an arbitrarily directed 
magnetic field, i.e., the coexistence of a transverse mag- 
netic field H x along the x-axis and a longitudinal one H z 
along the z-axis, the Euclidean action is written as|25j 
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where V is the volume of the particle, 7 the gyromag- 
netic ratio, x± = Xi(™2/mi) and If 11 the longitudinal 
anisotropy. The polar coordinate 9 and the azimuthal 
coordinate </>, which are the angular components of mi 
in the spherical coordinate system, can determine the di- 
rection of the Neel vector. Dot over a symbol denotes a 
derivative with respect to the Euclidean time r. 

The classical trajectory corresponding to the Euclidean 
action (1) is determined by the equations 
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1 A ll FIG. 1: The effective potential -V(6, <p = 0). 
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where n = m/(K^), x = Xx/(-^||7 2 )> K(z) = l H x(z), 
Vg = dV/ 89 and = dV/ d<f>. The inverted potential is 

V{9,4>) = K\\(-sm 2 9 + 2h x sin9cos<p + 2h z cos9 

+ sin2 4> + -7T cos2 cos2 
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where h x r x } = H x t z -\/H c and H c = 2Ku/m. In Fig.l the 
effective potential —V{9, <p = 0) is drawn. The minima 
of the potential correspond to the equilibrium orienta- 
tion of the Neel vector. The metastability condition that 
dV{9,Q)/d9 = and d 2 V(9,0)/d9 2 = determines the 
critical parameters at which the barrier vanishes. 

In the high-temperature regime the sphaleron solution 
of Eqs.|0 and is (#o,0o = 0)- is the position of 
the top of the potential barrier — V(6q,0), and is deter- 
mined by dV(9,0)/d9\g = and d 2 V(9,0)/d9 2 \g B > 0. 
9q has a cumbersome analytical form and its numerical 
result will be applied below to determine phase bound- 
ary lines. Furthermore, its behavior is illustrated in Fig. 2 
for given h x and y. Above, y = x/n 2 {= X±K\\/ rn ' 2 ) and 
the parameter y indicates the relative magnitude of the 
non-compensation. For large non-compensation (j/ < 1, 
i.e., m 3> y/x±K\\) and for small non-compensation 
(y 3> 1, i.e., m <C \/xaJ^\\)j the system becomes fer- 
romagnetic and nearly compensated antiferromagnetic, 
respectively. [13 Note that 9 = tt/2 for h z = 0. The 
crossover behavior of the escape rate of this model from 
quantum tunneling to thermal activation can be obtained 
from the deviation of the period of the periodic instan- 
ton from that of the sphaleron. To this end we expand 



FIG. 2: 9o versus h z for some given h x and y. 

(0,4>) about the sphaleron configurations 9q and cj> , i.e., 
9 = 9 + t)(t) and </> = <j) + £(t), where cj) — 0. Substi- 
tuting them into Eqs.© and one yields the following 
power series equations of the fluctuation fields 77 and £ 

( G{( V ,0 \,( Gl(r,,0 \ ( Gl( V ,0 \, _ Q 
I GJ(V,0 ) + { GUV,0 ) + { GS(r?,0 7 ' 

(5) 

where G?i, G2, G3, ...denote terms which contain linear, 
quadratic, cubic and higher powers of the small fluctua- 
tions, respectively: 
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It is introduced that Vgg = [d 2 V/d9 2 ]g = e ^ =( j, 

[d 2 V / d9d<j) 2 }e =eo ^=s , and so on. 

Denoting #Q(t) = (f?(r), £(t)), we have 50(t + f3H) = 
<5f2(T) at finite temperature and write it as the Fourier 
series SCI(t) = J2n=-oo SQ n exp[iLo n r], where to n = 
2im/ (3h. Since simple analysis shows that 77 is real and £ 
imaginary, to the lowest order we write them in the form 
77 ~ a0\ cos(<x>t) and £ ~ ia<p 1 sin(o;r). Here a serves as a 
perturbation parameter. Substituting them into Eq.JSJ) 
and neglecting terms of order higher than a, we obtain 
the relation 
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Next, let us write rj ~ a9\ cos(cjt) + ?y 2 , and £ ~ 
ia^x sin(wr) + i£ 2 , where ?7 2 and £ 2 are °f the order of 
a 2 . Inserting them into Eq.JSJ), we arrive at u> = c<j + and 

r/ 2 = a 2 po + cl 2 P2 cos(2a>r), £ 2 = a 2 q 2 sin(2u;T), (10) 



where the analytic forms of coefficients po, p 2 and q 2 are 
cumbersome, which are listed in the Appendix. 

This implies that there is no shift in the oscillation 
frequency. In order to find the change of the oscillation 
period, we proceed to the third order of perturbation 
theory by writing rj ~ a0i cos(wr) + i] 2 + r/ 3 , and £ ~ 
ia^i sin(wr) + it; 2 + i£ 3 , where r] 3 and £ 3 are of the order 
of a 3 . Substituting them again into Eq.lJSJ), and retaining 
only the terms up to 0(a ), we have 
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where 
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The forms oi gx{h x ,h z ,y) and g 2 {h x ,h Zl y) are 
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where w = nu + and A = <fi 1 /6x- Again, y = x/n 2 {= 
Xj_K\\/m 2 ) and the parameter y indicates the relative 
magnitude of the non-compensation. Also, pQ, p 2 and 
q 2 are obtained by replacing ip 1 by A0i and dropping 
2 in pq, p 2 and q 2 , respectively. It can be shown that 
for h z = Eq.lO is reduced to the case corresponding 
to uniaxial antiferromagnetic particle with a transverse 
magnetic field only has been investigated in Ref . |l7| . 

As shown by Chudnovsky 3] , if the oscillation pe- 
riod r is not a monotonic function of a, where a is 
a function of E in the absence of dissipation, the sys- 
tem exhibits a first-order transition. Thus, the period 
t(= 2-k/uj) in Ea. (|TT|) should be less than r + (= 2tt/oj + ), 
i.e., u> > lo + for the first-order transition. It implies 
that g(h x , h z ,y) > in Eq. Qllfl for the first-order tran- 
sition, and g(h x , h z , y) = determines the phase bound- 
ary between the first- and the second-order transition. In 
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FIG. 3: Phase diagram h x (y) for h z = 0,0.1 and 0.2. 



FIG. 5: Phase diagram h x (h z ) for y = 0.2 and 0.3. 



FIG. 4: Phase diagram h x (h z ) for y = and 0.1. 



this case the three parameters h x ,h Zl y should be treated 
simultaneously, which is not a simple problem. In the 
present work we will fix one parameter and then com- 
pute the boundary curve with the other two parameters. 
We first solve Ea. Ijllfl numerically to obtain the phase 
boundary lines h x (y) 's for several values of h z , which are 
plotted in Fig. 3. From Fig. 3, an immediate observation 
is that the first-order region for a given h z diminishes as 
y increases, which shows the same trend as the h z — 
case. Thus, it is evident that the region for the first-order 
transition is greatly reduced as the system becomes ferri- 
magnetic and there is no first-order transition in almost 
compensated antiferromagnetic particles. The result co- 
incides with Ref . [ri| . Fig. 3 also shows that with increas- 
ing h z the variety of first-order region is not a simple 
case. For the small y case with increasing h z the first- 
order region is shrunk, while for the larger y case the 
longitudinal field h z favors occurrence of the first-order 
transition. For instance, for the case of y = 0.05, the 
maximum values of the transverse filed h x for occurrence 
of the first-order transition are h x — 0.203, 0.195 and 
0.175 for h z — 0, 0.1 and 0.2, respectively. On the other 



hand, the first-order region vanishes beyond y ~ 0.46, 
0.55 and 0.89 for h z — 0, 0.1 and 0.2, respectively. 
This can be qualitatively understood from the consid- 
eration that the height of the effective potential barrier 
decreases as h z increases, whereas the height increases 
as y increases, therefore there is a competition between 
the longitudinal field and the relative magnitude of the 
non-compensation. When y — 0, the fact that the region 
for the first-order transition decreases as the longitudinal 
field increases results from a flattening of the peak of the 
barrier. [TH For the small y case (i.e., the large noncom- 
pensation), the crossover behavior of the ferrimagnetic 
system still keeps qualitatively that of the ferromagnetic 
one. However, for the larger y case the exchange interac- 
tion plays the role of effective magnetic field and so, for a 
given small transverse field the region for the first-order 
transition increases as the longitudinal field increases. To 
illustrate further the effect of the longitudinal field h z on 
quantum-classical transition, we next calculate the phase 
boundary lines h x (h z Ys for several values of y, which are 
shown in Figs. 4 and 5. In Fig. 4, for the case of y = cor- 
responding to uniaxial ferromagnetic system, the phase 
boundary line is plotted by a dotted line. Obviously the 
line coincides with Fig. 13 in Ref.jllj. in which quantum- 
classical transition in a uniaxial ferromagnetic system 
with a transverse magnetic field and a longitudinal one 
was investigated. For the case of y = 0.1, the phase 
boundary lines h x (h z Ys shift downwards and with h z in- 
creasing the critical value of h x decrease monotonically. 
Fig. 5 gives another case, in which the phase boundary 
lines h x (h z Ys show a kind of non-monotonic behavior. 
For instance, for the case of y = 0.3, the first-order re- 
gion vanishes beyond h x ~ 0.024 for h z — 0, while the 
maximum is h x ~ 0.112 for h z = 0.332. 

It was shown that quantum tunneling shall show up 
at higher temperatures and higher frequencies in antifer- 
romagnetic particles than in ferromagnetic particles of 
similar size. |25j Moreover, most ferromagnetic systems 
are ferrimagnetic, so nanometer-scale antiferromagnets 
are more interesting from experimental and theoretical 
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aspects. But a detailed comparison between the the- 
ory and experiment on quantum-classical transition re- 
mains a challenging task. It is very important to ob- 
tain the information on the magnitude of the quantity 
y for observing the first-order transition in real exper- 
iments. For the typical antiferromagnetic particle with 
X ~ 10~ 4 , Kh ~ 10 6 erg/cm 3 , and m ~ 500 emu/cm 3 , 
one can get the quantity y ~ 10~ 4 .pj| In this case, for the 
longitudinal field parameter h z — 0.4, the range of the 
transverse field parameter for observing the first-order 
transition is < h x < 0.127. It is note that Wensdorfer 
et al have performed the switching field measurements 
on individual ferrimagnetic and insulating BaFeCoTiO 
nano-particle containing about 10 5 -10 6 spins at very low 
temperature (0.1-6 K).[2q] Below 0.4 K, experimental 
results are quantitatively in agreement with the predic- 
tions of the macroscopic quantum tunneling theory with- 
out dissipation. The BaFeCoTiO nano-particles have a 
strong uniaxial magnetocrystalline anisotropy. 26J There- 
fore, the material is expected as a candidate to investi- 
gate quantum-classical transition of the escape rate of 
uniaxial ferrimagnetic or antiferromagnetic particles in 
an arbitrarily directed field. 

In conclusion, we have investigated quantum-classical 
escape rate transition for uniaxial antiferromagnetic 
particle with an arbitrarily directed magnetic field, 
i.e., the coexistence of the transverse and longitudinal 
magnetic fields. There are three parameters which can 
be controlled by experiment: relative magnitude of 
the non-compensation and two field parameters. The 
nonlinear perturbation method is used to obtain various 
phase diagrams for first- and second-order transition 
depending on the three parameters. It is shown that 
the effects of the applied longitudinal magnetic field on 
quantum-classical transition vary greatly for different 
relative magnitudes of the non-compensation. 

The work was supported by China Postdoctoral Sci- 
ence Foundation under Grant No. 2002032138, National 
Natural Science Foundation of China under Grant Nos. 
10174015 and 10234010, and Research Grant Council of 
Hong Kong. 

Appendix 

The coefficients in Eq.^UJ are deduced by using soft- 
ware Mathematica3.0. 
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